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Abstract 

When computed to one-loop order in resummed perturbation theory, the non- 
abelian Debye mass appears to be logarithmically sensitive to the magnetic scale 
g 2 T. More generally, we show that in higher orders power-like infrared divergences 
forbid the use of perturbation theory to calculate the corrections to Debye screening. 
A similar infrared problem occurs in the determination of the mass-shell for the 
scalar propagator in 2+1-dimensional scalar electrodynamics. In this context, we 
provide a non-perturbative approach which solves the infrared problems and allows 
for an accurate calculation of the scalar propagator in the vicinity of the mass-shell. 
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1 Introduction 



Significant progress has been achieved during the last few years toward the understanding 
of the infrared structure of high temperature QCDjl], pj, |3|. The prominent role of the soft 
energy scales gT and g 2 T has been recognized, and the collective nature of the dominant 
behaviour at the scale gT has been properly understood. (Here, g = g(T) is the coupling 
constant at the temperature T, and we assume that g <C 1 in the high-temperature, 
deconfined phase of QCD.) This led to a systematic description, in classical terms, of a 
variety of collective phenomena like screening, Landau damping, or color oscillations H. In 
the case of screening, it is known that, to leading order in g, the electrostatic interactions 
are screened, with a screening mass ~ gT, while the magnetostatic interactions are 
not screened 0, [5]. These properties are shared by abelian and non-abelian plasmas (see 
Refs. H |7|, H, |9f for the abelian case). 

Important differences occur between abelian and non-abelian gauge theories when 
corrections to the leading order Debye screening are considered. In the abelian case, 
perturbation theory can be used to calculate the corrections to the leading order Debye 
mass (see f|] and references therein). In QCD, infrared divergences occur in such a 
calculation, whose origin is the coupling of the chromoelectric field to the unscreened 



magnetostatic fields (see, e.g., Ref. |T0| for a survey of the computations prior to 1993, 
and also Refs. [11-15] for more recent calculations). For example, at one-loop order, there 
is a logarithmic singularity, widely discussed in the literature ]TI|, 0; H3- But we shall see 



that the difficulty is actually more serious, since power-like infrared divergences occur in 
the higher orders. 

The existence of infrared divergences invalidating the perturbative expansion of ther- 
mal QCD is well known for the magnetostatic sector, where power-like infrared divergences 



are indeed expected in higher order calculations of thermodynamical quantities jl6|, [L7|] . It 
has long been recognized that, because only static modes are involved, these divergences 
are essentially those of an effective three-dimensional theory |T8 . The divergences that 



we shall encounter here, which are also those of an effective three-dimensional theory, are 
of a slightly different nature. They occur in the perturbative evaluation of the polari- 
sation tensor of the electrostatic gluon on the tree-level mass-shell. Similar divergences 



are encountered in the calculation of the quasiparticle damping rates (see, e.g., |1£] and 
references therein). All such divergences could be removed by introducing an infrared 
cut-off A in the magnetostatic sector. However, this is not a satisfactory solution for at 
least two reasons. In QCD, there is a common belief that such a cut-off is indeed gener- 
ated dynamically in the form of a magnetic mass A ~ g 2 T |T7j, BO]. But for such a value 
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of A there are infinitely many terms in the perturbative expansion which contribute to 



the same order, a situation analogous to the Linde problem |T6j. The second reason is 
that we shall identify similar mass-shell divergences in the evaluation of the static scalar 
propagator in thermal scalar electrodynamics (SQED). And we know that there is no 
magnetic mass in abelian gauge theories || 

Thus, although we expect the picture of Debye screening to hold in higher order 
calculations, for reasons which will be detailed in the next section, it appears that the 
corresponding value of the screening mass cannot be computed in perturbation theory 
beyond the leading order. We are thus led to look for a non perturbative description which 
allows for the treatment of the large degeneracy of states involving massless magnetostatic 
fields. We shall propose such a treatment for SQED, and obtain the mass-shell behaviour 
of the scalar propagator without any infrared regulator. 

Our analysis relies on a non-perturbative approximation to the Dyson- Schwinger 
equations. The method that we use, known as the gauge technique, has been developed 
originally 0, ^] in relation to abelian gauge theories in four dimensions, and has been 
found to be particularly convenient for the study of the infrared structure of the prop- 
agator. It has the advantage to preserve the correct Ward identities, and the expected 
analytical properties of the propagator. Within this formalism, we shall be able to de- 
termine the infrared behaviour of the three dimensional scalar propagator. We shall find 
that the mass-shell singularity, which is a simple pole in leading order, turns into a branch 
point, whose location can be shown to be gauge-fixing independent. 

The plan for the rest of the paper is as follows. In section 2, we introduce the 
screening function, and discuss its analytic properties. In section 3, we critically analyze 
the previous computations of the Debye mass in the resummed one-loop approximation, 
and show that the loop expansion generates power-like infrared divergences. We identify a 
similar problem in the charged sector of SQED. In section 4, we present a non-perturbative 
approach which allows us to study the mass-shell behaviour of the scalar propagator in 
2+1-dimensional SQED. The last section summarizes the conclusions. 



2 The screening function 

In electrodynamics, the potential between two static charges in a medium can be calcu- 
lated from the electrostatic propagator in the medium, to be referred as the screening 
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function in what follows, 



ikx 



d^k e* 

S(x) = I dx D 00 (x , x) = / — k2 _ Uoo{ ^ k) , (2-1) 



where x = |x| is the distance between the two static charges, k — |k| = v k 2 , n o(0, k) is 
the static (k = 0) electric polarisation tensor, and -D 00 (0, /c) the electrostatic propagator: 
D 00 (0, k) = l/(k 2 — n o(0, k)) = S(k). We shall need later to analytically continue S(k) 
to complex values of k. In most cases to be discussed in this paper, S(k) will be obtained 
explicitly as an even function of k. When this is not so, we shall regard S as a function 
of ^/k 2 , i.e. S(k) = S(Vk 2 ), before doing the analytic continuation. 

A similar screening function describes the interaction of two static color charges in a 
quark-gluon plasma. In this case, however, the polarization tensor is gauge-dependent [p4] , 
E5L |TIfl . A fully gauge invariant treatment of chromoelectric screening should start from 
a gauge invariant object, such as the Polyakov loop. But, in perturbation theory, the 
leading long range behaviour of the correlator of two Polyakov loops is determined by the 
screening function (2.1) (see Refs. [T^, [T3|]). (We are not implying here that perturbation 
theory correctly describes the long range behaviour of the Polyakov loop correlator, which 
is presumably dominated by glue ball intermediate states f26j.) For this reason, we shall 
concentrate on this simpler object here. In fact, the long-distance behaviour of S(x) turns 
out to be gauge-independent. This may be expected from general arguments [0, and 
will be verified explicitly. 

Let us then return to eq. (|2.1|) which we rewrite as 



S( X ) = J- ^LkS{k)e ikx . (2.2) 
2irx J-oo 2ni 

In leading order, So(k) = l/{k 2 + m 2 ) where 

a 2 NT 2 

m 2 = -OT 00 (0,0) = ^— (2.3) 

is the leading-order screening mass squared and 511^ is the polarization tensor in the 
"hard thermal loop" approximation [1-5] (we consider here a pure gluonic plasma). The 
integral (|2.2| ) may be computed by continuing the integrand to complex values of k, and 
by closing the integration path in the upper half of the complex /c-plane. One then picks 
up the contribution of the pole k = im, and gets the familiar screened Coulomb potential, 
S ( x ) = e- mx /A-Kx. 

In abelian plasmas, the higher order corrections do not change significantly the 
picture. The singularity of S(k) which controls the asymptotic behaviour of S(x) remains 
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a pole on the imaginary axis; recently, this has been verified explicitly up to two orders 
beyond the hard thermal loop approximation ||. In such a case, we can define a screening 
mass mo as the solution of IB I 



m^ = -n 00 (0,A;) \ k , = _ ml . (2.4) 

This self-consistent equation admits also a meaningful perturbative solution^. 

Higher order corrections to Iloo(0, k) play a more dramatic role in QCD, where they 
alter the nature of the singularities of S(k). In order to discuss these corrections, it is 
convenient to remark that, in leading order, they can be considered as loop corrections 
in the effective three-dimensional theory obtained after integrating the non-static loops 
with static external lines (see [^SL [2^] and references therein). At the order of interest, 
and in Coulomb or covariant gauges, the corresponding Euclidean action reads 




where Di = di—igVTAi, = [D{, Dj]/ (igy/T), m is the leading-order electric mass from 
eq. ( |2.3| ), ( is the gauge fixing parameter, and 77 and fj are the ghost fields. The effective 
theory describes interacting static and long-wavelength (k ^ gT) fields. The magneto- 
static gauge fields A"(x) and the electrostatic field ^4q(x) are, up to normalizations, the 
zero-frequency components of the original gluonic fields. 

In the effective theory, Aq enters as a massive scalar field, whose propagator is the 
screening function. We shall write 

S-\k) = k 2 + m 2 + Z{k) (2.6) 

with £(&) denoting the self-energy corrections in the three-dimensional theory, i.e. we set 
n o(0, k) = —m 2 — E(&). It is easy to see that all the Feynman diagrams contributing 
to are analytic in k 2 for small momenta. Indeed, in any such diagram, one can 

choose the independent loop momenta so that the external momentum k flows only along 
the massive propagators. These can be expanded with respect to k, when \k\ <C m. In 
the resulting expression, the external momentum appears then only in the numerator, 
and rotational symmetry ensures that only the terms with even powers of k survive the 
angular integration. 

One may regard the effective action (|2.5| ) as the Euclidean version of a Minkovskian 
action in 2+1 dimensions. From this point of view, one expects S(k) to be analytic in the 
whole complex fc-plane, except on the imaginary axis. We shall assume that this property 
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indeed holds and write 



oo 



S(k) = i du -|^L . (2.7) 







The spectral density p{oj) may be calculated from the discontinuity of S{k) accross the 
imaginary axis: 

p (lo) = — lmS(k = i(u + ie)) (2.8) 

71 

with uj real. Because S is an even function of k, p(u) is an even function of u (this is 
easily verified by noticing that S(k) is real for k real, and applying the Schwartz reflexion 
principle: S*(k) = S(k*)). Thus, only the positive values of u are needed to represent 
S(k). At leading order, S(k) = l/(k 2 + m 2 ) and p(uS) = 2m5{uj 2 — m 2 ). We define 

E(u) = E(k — i(u + ie)) (2.9) 

where u is real. Then 

/ x 2u; ImS(w) . _ 

PH = - — t —\j T2- (2-10) 



7T 



uj 2 — m? — ReSfu;) 



Im T^iu) 



In a Minkovskian theory, one expects p{u) to be positive in a physical gauge. However, 
what is meant by a physical gauge is not the same in the present (2+l)-dimensional 
problem and in the original (3+l)-dimensional one. For the original problem at finite 
temperature, one can choose, as a physical gauge, the strict Coulomb gauge, that is, the 
limit £ —>■ in eq. ( |2.5| ). In the Minkovskian theory in 2+1 dimensions, this does not 
correspond to a Coulomb gauge, but rather to a Landau gauge which involves unphysical 
degrees of freedom. Since these latter do not give positive definite contributions to the 
spectral density, p(u) is not then necessarily positive for the theory defined by the effective 
action eq. (|2.5| ). 

Since the three-dimensional theory (|2.5|) is superrenormalisable, S(fc)//c 2 — > as 
k — > oo. It follows then from eq.(2.6) that, as \k\ — ► oo, S(k) ~ 1/k 2 . This property 
allows us to derive a sum rule for p. First we note that, owing to the asymptotic property 
of S just mentioned, we have 

t^ k S(k) = l, (2.11) 
J Ltxi 

where the contour is a circle at infinity in the complex A;-plane. Then we replace in this 
equation S(k) by its expression ( |2.7| ) in terms of the spectral function. The contributions 
of the two poles at k — ±ia; give 

/•DO 

duop(uo) = l. (2.12) 
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We shall see later that, in some approximations, S(k) may have poles away from the 
imaginary axis, which are not accounted for by the spectral function p(u). We shall argue 
later that such poles are unphysical, but it is nevertheless useful to keep track of them by 
writing 







The sum rule is then modified into 

/ dup(u) + y £(a i + aZ) = l. (2.14) 
Jo 

We conclude this section by summarizing the analytic properties that we expect 
the screening function S(k) to satisfy. We have given arguments suggesting that S(k) is 
analytic in the complex A;-plane, with singularities on the imaginary axis. Furthermore, 
in most gauges, S(k) is analytic in k 2 for small \k\, i.e. k <C gT. We shall verify that 
these properties are satisfied by the approximate S(k) that we shall obtain. We shall find 
that S(k) has branch cuts along the imaginary axis, starting at k — ±im* (see Fig. 1). 
The branch point dominates the asymptotic behaviour of S(x). According to eqs. ( |2.2| ) 
and (|2.7| ), we can write 

S(x) = — Hduj e- ux p(uj) , (2.15) 

A-KX Jm* 

which expresses the screening function as the Laplace transform of the spectral density. 
At large x, S(x) ~ f(x)e~ m * x , and we shall verify that m* is gauge-fixing independent. 

The previous analyticity arguments, which are sufficient to establish the exponential 
fall off of the screening function at large distances, may become invalid in some particular 
gauges. In particular, in the temporal axial gauge, the asymptotic fall-off of S(x) was 
found to be a power law rather than an exponential [|H|, |T5|j . It is likely however that this 



peculiar behaviour is a gauge artifact (see also Ref. [0). Indeed, the temporal axial gauge 



is known to lead to specific difficulties in the imaginary-time formalismjlT]]. It prevents in 
particular the power counting arguments leading to the effective three dimensional action 
T5D. 



3 Corrections to Debye screening are non-perturbative 

In perturbation theory, we expect the dominant singularity of S(k) to remain close to 
the leading-order pole at k = im. Thus, the determination of the Debye mass involves 



6 



the calculation of for k ~ im. Since E(A; = im) is infrared singular in perturbation 
theory, this leads to difficulties whose physical origin is analyzed in this section. We shall 
be led finally to the conclusion that perturbation theory cannot be used to estimate the 
corrections to the leading-order Debye mass. 



3.1 The polarisation tensor at next to leading order 

The one-loop graph contributing to £ is displayed in Fig. 2. It is readily evaluated as 

nk , m) ^ fNT j^J^^ D%M , (3 ,) 

where = qi/q) 

^•(0,q) = ^#^ + C^| (3.2) 
qz q^ 

is the free propagator for magnetic gluons. We are using here, and throughout, dimen- 
sional continuation in order to regularise the ultraviolet (mass) divergences. After com- 
puting the integrals, no UV singularity will actually subsist in the limit D — > 3. It is 
important to keep in mind that the limit D —>■ 3 will always be taken before discussing 
the infrared structure of the integrals. After a simple rearrangement, eq. ( |3.1|) takes the 
form 

2»™ f dD( l f 1 2 m 2 -k 2 

E(k; m) = g NT / - — ^\ — + — - — 3.3 

v ' ' y J (2tt) d \ q 2 + m 2 q 2 q + k 2 + m 2 v ; 



+ (C-l)(k 2 + m 2 



and gives, after an elementary integration, 



q • (q + 2k) 

q 4 ((q + k) 2 + m 2 ) 



E(fc;m) = am P ^arctan— + (C - 2) \, (3.4) 

I vixk m J 

where a = g 2 NT/47r has the dimension of a mass. 

The function ( |3.4|) has logarithmic branch points at k — ±im. The origin of this 
singularity may be seen on eq. ( |3.3| ): as k 2 — > — m 2 , the dominant contribution to the 
integral comes from the small q region, and when k = ±im the integral in fact diverges. 
To understand physically what happens, it is convenient to do a Wick rotation. We have, 
for real ui (recall eq. ( |2.9| )), 

^/ \ f [m 2 + uj 2 ) , m + us + ie ,, .1 . 

E(u;m ) = ami- ^ln + {(-2)}, (3.5) 

muj m — uj — ie \ 



7 



from which one gets 



~ CO 2 — TYl 2 

ImE(w;m) = na 9{oj 2 — m 2 ) . (3.6) 

co 



This imaginary part is proportional to 3>(w), the invariant phase-space for the decay of a 
particle of energy co into a particle of mass m < co and a massless particle. This is easily 
computed in the rest frame of the decaying particle as 



where e g = \/q' 2 + m 2 . Note that this 2-dimensional phase space does not vanish when 
u; — > m + , in contrast to the 3-dimensional one. This behaviour of the phase-space factor 
is responsible for the infrared divergences to be discussed further in section 3.3. 

Another noteworthy feature of eq. (|3.6|) is that ImE(w;m) > for co > 0, whereas 
with our conventions we would expect the opposite sign. The sign of Im £ is related to that 
of the spectral function according to eq. Q2.10I) , so that, in the one-loop approximation, 



p(u) is negative for all momenta co > m. In particular, from eq. ( |2.10| ) and the asymptotic 
form of the one-loop self-energy ( |3.5| ), namely £(u; — > oo) ~ iiraco, one obtains, for large 

CO, 

p{u) ~ -1— 2 . (3.8) 

CO* 

Thus, as alluded to after eq. (|2.10|) , in the present covariant gauge p(co) cannot be regarded 
clS db physical spectral density. 

Because the poles at k = ±im of the unperturbed propagator coincide with the 
branch points in the self-energy, which furthermore diverges in these points, the equation 
([2.4D cannot be used to calculate perturbatively the correction to the Debye mass. In 
fact the analytic structure of the propagator S(k) = l/(k 2 + m 2 + £(&)) is very different 
from that of the unperturbed one. It has branch points at k — ±im and, besides, a 
set of four simple poles at k — ±a ± ib, where the real numbers a and b are gauge- 
dependent To leading order in a, the values of a and b are given by a = a(n — 9), 
b - m (a/2)[C - 2 + ln(4m 2 /(a 2 + (b - m) 2 ))], with 9 = arctan(a/(6 - m)). It is 
instructive to follow the trajectory of these poles in the complex fc-plane, as a function 
of the dimensionless parameter a/m ~ g. In order to do so, we rewrite the inverse 
propagator as 

S-\k) = m 2 {l + x 2 + -f(x)\ (3.9) 
I m J 
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where x = k/m and f{x) = E/(am). For small coupling, the poles behave as indicated 
above. When the coupling increases, they follow the trajectories displayed in Fig. 3 (for 
the gauge Q — 2). There exists a critical coupling at which the poles become real. Beyond 
that, one of the pole flows toward m, the other being equal to an. We note that the latter 
regime corresponds to the small mass regime, which is attained here at strong coupling. 
The pole at k = an is the tachyonic pole already identified in the studies of the massless 
theory pO, p[. 



This analytic structure of the one-loop propagator, which contradicts the expected 
properties of S(k), leads to unphysical properties. Indeed, the relative magnitudes of b 
and m, which determine the asymptotic behaviour of S(x), depends on the gauge. If 
b > m, the long-range behaviour of the screening function remains dominated by the 
logarithmic singularity at z = im so that for x — > oo, S(x) ~ f(x)e~ mx /x. However, 
since the spectral function is negative, the pre-exponential factor f(x) is strictly negative 
(see eq. ( [2.15|) ), in contrast to the leading order result f Q = 1/An. If now b < m, the 
pole contributions dominate, and the screening function oscillates asymptotically. These 
changes of regimes for small changes in the parameters are physically not satisfactory. 



As we have mentioned after eq. ( [3.4|) , the singularities of the integral ( |3.3|) are de- 
termined by the small q region. They are therefore very sensitive to the small momentum 
behaviour of the magnetic gluon propagator. By allowing for a small gluon mass A, i.e. re- 
placing 1/q 2 — > l/(q 2 + A 2 ) in the integral (|3.3| ), one separates the mass-shell of the scalar 
particle and the threshold for gluon emission, and this removes the infrared divergences. 
A simple calculation gives then||ll|| (for A<m) 



„ x \2{m 2 -k 2 ) k m 2 + k 2 ) . 

E x (k; m )= ami -± — -Cretan — — - 1 + C - 1) , • 3.10 

[ mk m + X (m + X) z + k^ J 



In eq. ( |3.10| ), the branch point has now moved to m+X. In perturbation theory, S(k) has a 
pole at k — i(m + 5m) where 5m = H\(k = im; m)/2m ~ aln(2m/A). This perturbative 
analysis is consistent as long as the new pole does not move back into the cut, that is, as 
long as aln(2m/A) < A. However, if A ~ g 2 T, which is the order of magnitude expected 
for the magnetic mass, aln(2m/A) ~ g 2 Tln(l/g) 3> A, and one gets an inconsistency. 

One way to keep the pole separated from the branch cut is to change m in eq. (|3.10| ) 
into m£>. That puts the branch point at m^ + A, whatever the value of is. The Debye 
mass mn is then obtained by solving self-consistently the equation 

m 2 D = m 2 + T,\(k = imp ; mo) ■ (3-11) 

The pole at k = imn remains below the branch point at k = Ojud + X), so that it controls 
the long range behaviour of the screening function. However the sign problem alluded to 
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earlier is not solved. The sign of S(x) at large distances is determined by the residue at 
the pole. A simple calculation gives 

2\ e ~m D x 

S{X) ~— 2A + (C-3)c ^' (3 ' 12) 

the first fraction being the residue just mentioned. Thus, in Feynman's gauge for instance, 
S(x) becomes negative if A < a. Note that, in perturbation theory, one would expect 
the residue to be close to unity; the above formula shows that this only happens if the 
infrared cut-off is large enough, i.e. A 3> a, which is not to be expected. The above 
procedure leading to eq. Q3.11| ) is an attempt to go beyond perturbation theory, which 



ignores, however, all the vertex corrections. The fact that the latter may be important 
is suggested by the non-perturbative character of the residue. Nevertheless, we shall see 
in section 4.4 that, in the presence of a magnetic mass, the result ( 3.11J) remains correct 
even when vertex corrections are taken into account. 

In closing this subsection, let us mention that, for nonvanishing A, the same correc- 
tion to the Debye mass as obtained above, i.e. 5m w a ln(2m/A), can be deduced from the 
long range behaviour of the correlator of two Polyakov loopsfl^, Note however that 
the two calculations are not independent since, in perturbation theory, the asymptotic 
behaviours of the Polyakov loop correlator and of the screening function involve the same 
integrals. 



3.2 A similar problem in scalar QED 

In the high temperature limit, and at leading order in e, the static and long- wavelength 
(k < eT) correlation functions of scalar QED can be calculated from the effective three- 
dimensional euclidean action || 

Seff = J d 3 x (~ F* + ^{d^f + 1 (Mo) 2 + \ ml Al 

+(A0) t (A0) + m 2 ^(f) + +e 2 TA 2 ^(f) + ^(j)) 2 ) . (3.13) 



4 

The new notations here are as follows: 0(x) is the complex scalar field, m 2 el = e 2 T 2 /3 is 
the leading order electric mass, and m 2 = e 2 T 2 /4 is the charged particle thermal mass. 
We are only interested here in the mass-shell singularities associated with the interaction 
between the charged particles and the massless transverse photons. We shall therefore 
restrict ourselves to the sector of (|3.13j ) which describes these interactions: 



3 e ff = J d 3 x Q F 2 + ^-(diAi) 2 + (A0) f (A0) + m 2 0V) • (3.14) 
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There is an obvious similarity between this action and the corresponding one for the hot 
QCD plasma, eq. ( |2.5|) : they both describe massive charged particles (Aq in QCD and 
(f) in SQED) in interaction with massless gauge fields. In particular, to one loop order, 
the scalar self-energy E(fc;m) is also given by eq. ( |3.4| ) ||, and most of the discussion in 
section 3.1 applies to SQED as well. Of course, essential differences persist between the 
two theories in the dynamics of the gauge fields themselves. In particular, in SQED, the 
transverse photons remain massless to all orders in perturbation theory || |J, so that one 
cannot invoke anymore a magnetic mass to regularise the mass-shell singularities, as we 
did in eq. (|3~ToD for QCD. 

The massless version of the theory ( |3.14| ) has been studied extensively [IE]. In this 



case the one-loop scalar self-energy generates a tachyonic pole at k ~ a. This can be 
seen from the expression (|3.4j ): in the limit m/k — > 0, = — irak, and the inverse 

propagator S(k) = k 2 — irak vanishes at k = na. The infrared behaviour is improved by 
the resummation of the one-loop polarization tensor in the internal photon line in Fig. 2. 
To see that, consider the photon polarization tensor 11^(0, q) = (5y — g^g J )n r (0, q). At 
one- loop order in the effective theory we have[]3~3|, |9| 



I * M = B "(^" , *"s;- 1 l- (315) 

In the massless limit m/q — > 0, ITj^O, q) = na q/A is linear m. q, so that it dominates over 
the contribution of the bare inverse propagator at small momentum, that is, as q — > 0, 
Dij(0,q) ocl/q. 

Although this softening of the photon propagator does not solve entirely the in- 
frared problem, it is clear that the loop insertions in the internal photon lines do play an 
important role in the massless case. This is not so in the massive theory. The reason is 
that, when m^fl, and to all orders in perturbation theory, the polarization operator is 
expected to vanish at least as q 2 when q —>■ 0|| (in particular, Ely (0, q) ~ (a/6m)q 2 
as q — > 0); thus, in the massive theory, the low momentum behaviour of the resummed 
photon propagator is not different from that of the bare propagator. 



3.3 The need for a non perturbative treatment 

The infrared divergences that arise in the one-loop calculation signal, in fact, a breakdown 
of perturbation theory. There exists indeed an infinite number of multi-loop diagrams con- 
tributing to which become infrared singular as the external momentum approaches 
the tree-level mass-shell, i.e. when k 2 — > —m 2 . The one-loop diagram represented in 
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Fig. 2 is logarithmically divergent as k — > im. Consider the two loop diagram of Fig. 4a. 
Its complete infrared behaviour is calculated explicitly in the Appendix, but the leading 
terms can be obtained by simple power counting. The most divergent contribution is in 
the integral 

d 3 qd 3 p 

(3.16) 



(g 2 + A 2 )(p 2 + A 2 )(k-q) 2 (k- (p + q)) 

where we have added a small mass to the photon in order to facilitate the power counting. 
The integral over q is linearly divergent as A — > 0. The same result holds for the diagram 
in Fig. 4b, which involves vertex corrections, and it can be verified that the leading di- 
vergences of the diagrams 4a and 4b do not mutually cancel. A similar power counting 
argument can be extended to all Feynman diagrams involving no correction to the mag- 
netic photon line, such as the one displayed in Fig. 5. The result of power counting is 
that, close to the mass shell, a n-loop graph (n > 2) diverges like (a/A) n_1 , up to powers 
of ln(a/A). 

Physically, the origin of the infrared divergences is the degeneracy between the 
mass-shell of the charged particle and the threshold for the emission of n (n > 1) massless 
transverse photons. Then, the determination of the mass shell requires solving the theory 
in the subspace of these degenerate states. Naively, one would expect the coupling to two 
or more photons — which brings in more powers of g — to be less important than the 
coupling to a single photon. This is what happens in 3+1 dimensional electrodynamics. 
In the present case, the low dimensionality of the space-time amplifies the effects of the 
degeneracy, in such a way that the couplings to any number of photons become equally 
important. 

Similar divergences arise in QCD as well. Some of the relevant diagrams are actually 
the same as in SQED (e.g., Figs. 4a,b and Fig. 5), the scalar line in these diagrams being 
interpreted as an electrostatic gluon. Besides, there exist new divergent graphs involving 
the self- interactions of the magnetic gluons (see Fig. 4c for an example). The same power 
counting as above leads again to the conclusion that n-loop diagram diverge as (a/A)" -1 . 
One may argue that the infrared divergences are cured by the dynamical generation of 
a magnetic mass A. However, for A ~ g 2 T ~ a, as commonly expected [{Hq 1 , 17], all the 
aforementioned diagrams contribute to the same order in g. 

To summarize, the analysis of this section suggests that non perturbative methods 
are necessary in order to determine the correct mass shell behaviour. Such a method will 
be presented in the next section for the case of SQED. 
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4 An integral equation for the spectral density 



We present now an approximate, but non-perturbative, solution of the Dyson- Schwinger 
equation of scalar electrodynamics, which provides the behaviour of the scalar propagator 
near the mass shell. To this aim, we establish a linear integral equation for the spectral 
density p{ui) using the so-called gauge technique^^. This equation performs a partial 
resummation of the most infrared singular diagrams in a gauge-invariant way, i.e. by 
respecting the Ward identities. When applied to four- dimensional abelian gauge theories, 
it provides the correct mass-shell behaviour for charged particles p2[ . 



4.1 The quenched approximation 

The four skeleton diagrams which enter the Dyson- Schwinger equation for the scalar 
propagator are displayed in Fig. 6. It can be verified by power counting that, at a given 
order in e, the most singular diagrams are obtained from the perturbative expansion of the 
first graph, Fig. 6. a, where we keep the photon propagator at the tree-level (see Fig. 7). 
These are precisely the diagrams discussed in section 3.3. Thus, the Dyson-Schwinger 
equation that we wish to solve, and to which we refer as the quenched approximation, is 

E(k) = -e 2 T J ^(2* i + ft )£>J(q)r i (k+q ) k)5(k + q). (4.1) 
In this equation, S and Tj are the full propagator and vertex, related by the Ward identity: 

q t r,(k, k + q) = S-\k + q) - S'^k) . (4.2) 

The most general vertex function which is consistent with this identity and which is free 
of kinematical singularities is of the form[[3f| 

r 4 (k, k + q) = (2 ^ qf {S-\k + q) - S^k)) + A(6 tJ - q^kj , (4.3) 

where A = A(k 2 , q 2 , k ■ q) is an unknown scalar function. According to the usual termi- 
nology in the literature, we shall refer to the two terms in the r.h.s. as the longitudinal 
and transverse pieces of the vertex function, respectively. The second term, proportional 
to A, is indeed transverse to the photon momentum q. However, the first term is not 
parallel to q: it involves a non-trivial transverse piece which is completely determined by 
the Ward identity (see also eq. ( f4.6| ) below). At leading order, S'~ 1 (k) = k 2 + m 2 , A = 0, 



and Tj(k, k + q) given by eq. ( |4.3|) reduces to the tree- level vertex, 2fcj + q^ 
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As q — > 0, eq. (O) must reproduce the differential form of the Ward identity, 
namely 

Together with eq. (|4.3| ), this implies that A — ► as g — > 0, for any /c. This suggests that 
only the longitudinal vertex becomes singular as q — > and /c 2 — ► —m 2 D ] the corresponding 
singularities are, of course, those of the 2-point function. This conclusion is supported by 
perturbative calculations^]. Thus, in order to get the leading mass-shell behaviour, we 
need only consider the longitudinal vertex. That is, we make the following ansatz for Tj : 

r,(k, k + q) = J h + * (S-\k + q) - S-\k)) . (4.5) 



[2k + q) ■ q 

As anticipated, this has a component transverse to q: 

(<% - g,g,)r,(k, k + q) = 2( £ k ~™^ {S-\k + q) - S-\k)) . (4.6) 

It is this component which couples to the physical, transverse piece of the photon propa- 
gator (|3.2|). The longitudinal piece of (|4.5|) couples only to the gauge degrees of freedom 



ofS&(q). 

With the ansatz (|4.5| ) for Tj , the Dyson-Schwinger equation (|4. 1]) becomes 

= - e2r / <& B «< q) (2 *'(2k+q)'r j> (i - s " (k)s<k + q> ) ■ <47) 

We are eventually interested in the self-energy close to the mass-shell, £(& — > imp). By 
definition, S~ 1 (k = imp) = 0, so that the term proportional to S' _1 (/c) in the r.h.s. of 
eq. (4/7) does not contribute to £(& = imp), unless the integral over g diverges. Since this 
integral is indeed potentially divergent, we introduce a regulator, in the form of a small 
photon mass. In doing so, it is important for what follows to keep explicit the distinction 
between the physical and the unphysical states in the photon propagator. The general 
structure of the exact magnetostatic propagator in finite-temperature four- dimensional 
SQED is§ |: 

with H"r(0, q) = IIjj(0, q)/2. We have already mentioned, in section 3.2, that n r (0, q) oc q 2 
as q — > 0. However, as a convenient IR regularisation, we give the transverse photons a 
mass and set temporarly IIt(0, q) = \m (\m "C rn). We also regularise the gauge sector 
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of Dij by replacing 1/q 2 — > l/(q 2 + A 2 ) in the second term of eq. ( |4.8| ). Both ( and 
A should disappear in the evaluation of physical quantities. Then we get from eq. (|4.7|) 
E = E^ + £<5E where 

Ei( *> = - e2T IS^TTW, (^T& (' - s " (k)s(k + q » • (4 ' 9) 

is the self-energy in the Landau gauge, while 

is the gauge dependent part of E. Note that the above equation for <5E is actually 
independent of the ansatz used for the vertex Tj, since it follows directly from eq. (|4.1| ) 
and the Ward identity (PO). 



Consider now eq. (|4.9| ) in the on-shell limit k — > zmx>. As long as we keep the 
infrared regulator A at 7^ 0, the integral is convergent and the term proportional to S^ 1 (k) 
vanishes on the mass-shell. Thus 

d D q 1 4(k-q) 2 



E(& = im D ) = —e 2 T J 



(2tt)^ q 2 + \ 2 M (2k + q)-q 



, 2m n 
~ 2am D \n— — . (4.11) 

k=imo 



On the other hand, in the physical limit Am — > 0, not only does the estimate (|4.11 



become logarithmically divergent, but the integral multiplying S 1 (k) also diverges on 
mass-shell, since, in this limit, it is proportional to 

f ^ (4 12) 

J q 2 (q 2 + 2k • q)(g 2 + 2k • q + E(k + q) - E(Jfe)) ' 1 ' ' 

Thus the use of an infrared regulator does not allow us to explore in a simple way the 
behaviour of the scalar propagator near the mass shell. For this, more powerful technics 
are needed, such as that developed in the next subsection. 

Consider finally the gauge-dependent part of the self-energy, i.e., <5E, eq. (|4.10|) . 

Since 

d D q q-(2k + q) _ f d D q q 2 3A 



(2n) D (q 2 + A 2 ) 2 J (2tt) d (q 2 + A 2 ) 2 8tt ' 
which vanishes as A — > 0, the r.h.s. of eq. ( |4.10|) reduces to 

ffi(t) = «T^(k) / g D ^ ff g(k + q). (4.14) 

The usefulness of A appears when considering the on-shell limit of this equation. If we 
set A = 0, then the g-integral diverges on the mass-shell, and the limit 5E(/c — > irrin) 



15 



is not obvious. But if we keep A ^ 0, we obtain <5£(fc = imp) = 0. This guarantees 
the gauge-independence of the mass-shell, since the equation S~ 1 (k = im D ) = reduces 



then to m 2 T 



D 



m 



+ Sl(/c 2 = —m 2 D ). This procedure for taking the on-shell limit in 



the presence of an infrared regulator has been proposed in Ref. 
resummed one-loop approximation. (See also Refs. 
the computation of the damping rates.) 



11], at the level of the 



33| for a similar problem in 



4.2 The integral equation 



We shall transform now the Dyson-Schwinger equation (|4.7| ) into an integral equation for 
the spectral density p{uj). In this way, we take automatically into account the expected 
analytical properties of both the propagator and the vertex. First we put eq. ( f4.1|) in the 
form 

1 = (k 2 + m 2 )S(k) - e 2 T J -^L (2 kl + ft )Z>° (q) S(k + q) T,(k + q, k) S(k). (4.15) 

Then we use the spectral representation Q2.7|) of the propagator to rewrite the ansatz (|4~5l ) 
for the vertex function as 



S{k) r 4 (k,k + q),S(k + q) 



(2kj + gj)p{uj) 
n (uj 2 + k 2 )(u; 2 



k + q) 



(4.16) 



When eqs. (|2.7p and (|4.16j ) are used in ([4.15|) , the following equation for p is obtained]^! 



+ k 2 



k 2 + m 2 + £(fc; s) 



(4.17) 



where S(/c; s) is the one-loop expression ( |3.4| ) in which the mass m is replaced by s. That 
is. 



isk 



s — ik 



(4.18) 



where a = e 2 T /Aw. 



As it stands, eq.( f4.17| ) is not easy to solve. To make progress, we do a Wick rotation 
k ^ i(u + ie) to time like momenta. The integral equation becomes then 



d.s 



s 2 — uj 2 — it 



-u) + m + Yi{&\ s) 



(4.19) 
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By taking the imaginary part of ( [4.17 ) we obtain an equation which is linear in p and 
homogeneous: 

— p(u) \oo 2 -m 2 -Ret(uj;uj)} = f°° ds - lmt(uj;s) . (4.20) 

2uj 1 1 Jo s 2 — uj 2 

At this point, all what we have done applies to 3+1 dimensional SQED as well. 
It is interesting to see what happens in this case to enlighten the difference with the 
2+1-dimensional case. In 3+1 dimensions, 



ImS( W ; s) = ^ J )e2 ^±^- (uj 2 - s 2 ) 6(uj 2 - s 2 ). (4.21) 

l07T UJ Z 



The real part, Re£(u;;u;), is regular after UV renormalisation and it can be combined 
with the parameter m 2 in the l.h.s. of eq. (|4.20|) to define the physical mass (which we 
continue to denote by the symbol m for simplicity). The equation for the spectral density 



is then [[22 



(to — m )p(uj) 



(C-3)e 2 cu 

r2 



£ ds p (s )(l + il). (4.22) 



For u ~ m, this equation can be transformed into a differential equation which is easily 
solved. To do that, set F(w) = {u 2 — m 2 )p(uj), and verify that when u — > m, F(u) 
satisfies 

dF (C-3)e 2 2mF(u) 



du 87r 2 uj 2 — m? 

The solution of this equation gives the behaviour of p(uS) for uj —>■ m: 



(4.23) 



, MK _i_(^)^. (4.24) 
This is the correct behaviour, as obtained by a variety of other methods (see [031 and 



Refs. therein). 

An essential step in the previous calculation is the cancelation of the singularity at 
s 2 = uj 2 of the integrand in the r.h.s. of eq. ( [4.20 ) by the phase space factor {uj 2 — s 2 )/uj 2 



contained in ImEfu; s), eq. Q4.21 ). This phase space factor is very much dependent upon 



the dimension. Recall that, in 2+1 dimensions, it is simply ~ l/u> (see eq. (|3.7|)), so that 
in this case the imaginary part of the one-loop self-energy does not vanish at threshold; 
we have indeed 

2 2 

Im t(uj; s) = ixa - +S 9(uj - s) . (4.25) 

UJ 



17 



It follows that the integrand in the r.h.s. of eq. ( |4.20| ) is singular as s — > uj, suggesting 
that that eq. (|4.20|) is ill defined as written. However, according to eq. ( |3.5|) , Re£(c<j; s) 
is also logarithmically divergent when evaluated for s — > u, and it turns out that this 
divergence precisely cancels that of the r.h.s. To see this, we write the real part of the 
one-loop self-energy as 

171 UJ 2 + S 2 

ReS(cj;m) = (am + 2a — / ds — -. (4.26) 

uj Jo m 2 — s z 

The last integral diverges at its upper limit when m — > uj. When it is added to the integral 
in the r.h.s. of eq. (|4.20|) , the limit m —>■ uj becomes well-defined. The equation for p 
takes then the form 

2 2 

p(u) (uj 2 - m 2 - C a u) = 2a jT ds (p(u) - p(sj) . (4.27) 

This is our main equation. It has a number of properties which are worth emphasizing. 

i) The integral over uj in the r.h.s. now extend over a limited range of values of uj. This 
makes it well suited to study the behaviour of p(u) near threshold. 

ii) It is a linear integral equation for the spectral density. Being also homogeneous, it 
determines p only up to a constant factor, which in principle is fixed by the inhomogeneous 
equations (|4.17|) or (|4.19|) . It can be verified in particular that any normalisable solution 



of the integral equation Q4.17D satisfies the sum rule 

dujp(uj) = 1. (4.28) 



o 



iii) In the limit a — > 0, eq. ( [4.27 ) admits as a solution the free spectral function, p(uj) = 
2m5(uj 2 — m 2 ). 

iv) The asymptotic form of p(uj) for uj ^> m has been obtained from perturbation theory 
(recall eq. flOD ): 

a 

p(uj)~-2— . (4.29) 

UJ Z 

This is consistent with the integral equation ( [4.271 ) . Indeed, by assuming that the solution 
p(uj) falls off rapidly enough in order to be normalisable, one can obtain from eq. (|4.27| ) 
the following relation: 



uj 2 Jo 



a 

p(uj)~-2— I dsp(s). (4.30) 



When combined with the sum rule Q4.28j ), this equation yields the asymptotic behaviour 



Q4.29|) . It also shows that p must change sign. 
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4.3 Solving the integral equation near the mass shell 



From the analyticity of the Feynman diagrams (see the discussion after eq. ( |2.6D ), one 
expects p(u) to vanish identically in the vicinity of uj = 0. We show now that this is 
indeed the case for the solution of eqs. flilTD and fl£27D . Note first that eq. (|4~T7|) is valid 
for any k. Setting k = 0, and using E(& = 0; uj) = (au>, we obtain 

1 = r d uj^4- {™ 2 + (auj). (4.31) 

JO UJ ^ ' 

For this integral to converge, p(uj) must vanish sufficiently rapidly when uj —>■ 0. Then, 
assuming p(uj) to be regular, if it is not zero it is either increasing or decreasing for uj 
small enough. Assume, for example, that p is increasing, so that it is positive for small 
uj. Then, the l.h.s. of the integral equation ( 4.27 ), ~ —m 2 p(uj), is negative, while the 
r.h.s. is positive. One runs into a similar contradiction if one assumes instead that p is 
decreasing for small uj. The only acceptable possibility is that there exists m* > such 
that p(uj < m*) = 0. It may be furthermore verified on eq. (|4.27|) that p(uj) cannot have 
an isolated, 5-type singularity at uj = m*; that is, the spectral density is non- vanishing in 
any upper vicinity of m*. 



Because p(uj) vanishes when uj < m*, one can expand the integrand in eq. (|4.17j ) 
for small k without generating infrared singularities. In this way, one obtains sum rules 
involving higher and higher moments of 1/uj. For example, in the Landau gauge ( = 0, 
we have 

P{U) (4.32) 

(4.33) 

UJ I 

Such sum rules suggest that p is positive when uj 

To study the behaviour of p{uj) for uj ^ m*, we divide the s-integration in eq. ( |4.27| ) 
in two parts: from to m*, where p(s) = 0, and from m* to uj. After a simple calculation, 
the integral equation is rewritten as 

( 2 2 u + m * 

p{uj) [uj — m —C,auj + 2am —2a uj In 




uj — m' 



= 2a / ds^—- 2 (p(uj)-p(s)). (4.34) 

Jm* UJ — S x ' 

Assume that p(uj) is positive near the threshold, in conformity with the sum rules above. 
As uj — >• m*, the l.h.s. is dominated by the singularity of the logarithmic term, 

2m* 

l.h.s. w -2am* p(u) In , (4.35) 

uj — m* 
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which is negative (and gauge-independent). The r.h.s. must be negative as well, and this 
requires p{u) to be decreasing. Since p is positive and decreasing, it cannot vanish at 
m*\ p(m*) > 0. In fact, p is divergent at threshold, for, if p(m*) were finite, the integral 
in the r.h.s. of eq. ( |4.34| ) would vanish as uj — > m*, while the l.h.s., eq. ( |4.35| ), would be 



divergent. It follows that, close to the threshold, the integral in the r.h.s. of eq. ( |4.34| ) is 
dominated by the singularity of p(s) as s — > m*, so that we may approximate 

r.h.s. « — 2a / dsp(s). (4.36) 



u — m 

In writing this equation, we have neglected p(u) (which is finite as long as u > m*) next 
to p(s), and we have made the appropriate replacements, e.g., u + m* — > 2m*. From 
eqs. (|4.35| ) and ( 4.36|) , we obtain the following approximate form of the integral equation: 

1 r x 

xp(x)\n- = / dyp(y), (4.37) 



x Jo 

with the notation x = (u — m*)/2m* and the obvious identification p(x) = p(u = m* + 
2m* x). Note that a and all the other parameters (namely, m and Q have dropped from 
this equation. Its solution is easily found as 

p{x) = 2^ xlLy ' (4 - 38) 

where Z is a dimensionless constant. As expected, this is divergent as x —>■ + , but the 
divergence is integrable, as required for a solution of the integral equation. In the original 
variable, 

Z6(uj-m*) . 
p(u) « 2> uj^m*. (4.39) 

(«-m*)(ln ^) 

This is our main result. 

Near the mass-shell, the corresponding Minkowski propagator is (oj > 0) 



This is qualitatively different from the one-loop result discussed in Sect. 3.1. The in- 
verse propagator S^ 1 ^) vanishes at the mass shell, but its derivative is divergent there, 
a situation which is reminiscent of that in four dimensions ^3|. Thus, the integral equa- 
tion ( |4.15| ) provides, in (2+l)-dimensional SQED, a mass-shell behaviour for the charged 
particle propagator which is analogous to that obtained at the one loop level in (3+1)- 
dimensional SQED. In particular, the self-energy corresponding to the propagator (|4.40| ) 
behaves near the mass-shell as 

EMoc(^-m* 2 )ln m *~^~* e . (4.41) 
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It contains a factor uo 2 — to* 2 (compare, in this respect, with eqs. ( |3.6| ) and ( 4.21| )) which 



makes the imaginary part of ( |4.41|) vanish and change sign at the mass-shell. 

According to eq. (|2.15|) , the long-range behaviour of the screening function is de- 
termined by the spectral density near threshold. With p(uo) given by eq. (|4.39|) , we get 
S(x) ~ f(x)e~ m * x , so that to* plays the role of the screening mass. Note that if the 
integral equation (|4.37|) specifies the correct mass-shell behaviour, it leaves to* arbitrary. 
In principle, the value of m* could be obtained by solving the integral equation ( f4.17p for 
p(uj) for all the values of uo. But the approximations underlying this equation are only 
valid at small momenta. 

Now, if we cannot specify the value of to* anyfurther, we can verify that it is 
independent of the choice of the gauge parameter. To this aim, consider the gauge- 
dependent contribution to the scalar self-energy, as determined by eq. ( f4.14|) . This can 
be written in terms of the spectral density as 

Sm = e*T S-(k) / d.p M / ^ ^0 ^ +i l + qr , (4.42) 

and is readily evaluated, with the result 

fE W = .^(k)/«L, ^L_ , (4.43) 

If we were to set A = 0, we would obtain a non-vanishing contribution on the mass-shell: 
5E(k) = a S-\k) J duo — « rn* . (4.44) 

When taking the limit k — > im* in the above equation, we use the fact that, close to 
the mass-shell, the integral is dominated by the singularity of the spectral density at 
uo = m*. In order to get rid of the unwanted contribution ( [4.44| ), which would make the 
mass-shell gauge-dependent, we are thus led to keep A ^ before taking the on-shell limit 
(see the discussion after eq. (^4.14|) ). With this procedure, the coefficient Z in eq. (|4.39| ) 



is gauge dependent, and also A-dependent, as we explain now. Consider the equation 
S~ l = S^ 1 + C<5£ (with Sl the propagator in the Landau gauge) in the vicinity of the 
mass-shell, where the approximate form fl4.40|) holds (with Zl replacing Z in the case of 



Sl). By using eq. ( [4.43| ) for 5S(fc — > im*), one derives 

1 1 C" f°° , uop(uo) , A A . 

As A — > 0, the integral in the r.h.s. is essentially the on-shell propagator, so that it is 
divergent. Thus, for A small enough, the integral is dominated by the singularity of p(uo) 
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as uj — > m*, where the approximation ( |4.39j ) can be used. One deduces that, as A — > 0, 

(4.46) 



eq. ( 4.45|) takes the form 

11 



Z Z L 2Aln(2m*/A)' 
This equation determines the dependence of Z upon ( and A, in the limit A —>■ 0. We 
recall here that Zl is trivially independent of ( and A, since it is determined by eq. ( 4.17| ) 
with C = 0. According to sum rules like fl4.32|) -( [4\33| ), we expect Z L to be positive. 



4.4 The case of a non-vanishing magnetic mass 

For the sake of comparaison with previous computations, in particular those presented in 
section 3.1, it is instructive to consider the integral equation for the spectral density in 
the presence of an infrared regulator )i€m. Since the question of gauge invariance is 
not an issue here, we use a single regulator, in contrast to what we did earlier in section 
4.1. 

With A 7^ 0, the real and imaginary parts of £a(u;; s) are obtained from eq. ( |3.1C| ): 

Re Sa^; Co 1 ) = 2auj In — aiu , 

A 



uj 2 + s 2 „/ . ,.\ „ , na s(uj 2 — s 2 ^ 



ImS A ( W ; S ) = 7 ra^——d(uj-(s + \))-((-l)^'^ — S(u — (s + X)) 



(4.47) 



uj \ ' 2 uj 

Then, the homogeneous integral equation takes the form 

/ 2uj\ ru—X ijj 1 _)_ _g 2 

p(uj) I uj 2 — m 2 + auj — 2aa;ln — ) = (( — l)aujp(uj — A) — 2a / ds — -p(s) . 

V A / Jo uj — s 

(4.48) 

Note the following differences with respect to the case A = 0: i) ReSA(^;^) is finite; ii) 
the s-integration in eq. (|4.48 ) is now restricted to s < uj — A, thus avoiding the singularity 



of the integrand at s = uj; in) the gauge term is proportional to p(uj — A) and reflects the 
spurious pole at k — i(m + A) arising in the gauge piece of the self-energy ( |3.10| ). 

From eq. (|4.48j ), it is easy to establish that the mass-shell corresponds to a 5-type 



singularity, that is, to a simple pole in the corresponding propagator. To see this, consider 
eq. ( [4.48 ) for some uj satisfying m* < uj < m* + A, where m* is the mass-shell position; 



then, the r.h.s. of the integral equation vanishes since p(s) = 0iors<uj — \ < m*. The 
equation becomes, in this momentum range, 

p(uj) (uj 2 -m 2 + Ret x (uj;uj)) = 0, (4.49) 
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and it is solved by p(u) oc 8{u — m*), with m* determined by 



m * 2 = m 2 + ReE A (m*;m*). (4.50) 



Since Re £a(^; u), eq. Q4.47|) , is gauge independent, so is m*, the solution of the above 



equation. Eq. ( |4.50|) is identical to eq. (|3.11|) , obtained after a partial resummation of the 
scalar propagator which amounts to the replacement of the leading order mass m by the 
exact mass m D , that is, by using the propagator 

S(k) ~ - - 1 2 , (4.51) 
k 2 + m 2 D 

in the calculation of E(/c). The fact that vertex corrections do not seem to play any role 
in the determination of the mass shell is, strictly speaking, illusory. In fact, the true 
propagator in the vicinity of the pole is not ( [4.51| ), but rather z/(k 2 + m 2 D ). The residue 
z may differ significantly from unity; in principle, it could be determined if we were able 
to solve the full integral equation ( [4.4£| ). However, because of the Ward identity, the 
residue enters also the vertex correction, in such a way that it cancels against that of 
the propagator when the self-energy is computed at the mass-shell. It is therefore not 
needed to determine m*, but it enters as a preexponential factor in the asymptotic form 



of the screening function. We see now that the preexponential factor in eq. ( |3.12| ) is not 
consistently determined. 



5 Conclusions 

We have shown that the corrections to the Debye mass in high temperature non abelian 
gauge theories can be analyzed as corrections to the mass-shell of particles coupled to 
massless magnetic modes in a 2+1 dimensional effective theory. If one attempts to calcu- 
late perturbatively the Debye mass within this effective theory, one encounters power-like 
infrared divergences which signal a breakdown of the perturbative expansion. The possi- 
ble existence of a magnetic mass in QCD does not remove the essentially non perturbative 
character of the corrections. 

We have shown that similar mass-shell singularities occur in the evaluation of the 
scalar propagator in scalar electrodynamics. For this case, we have presented a non 
perturbative approach which allows for a complete description of the mass-shell behaviour. 
The mass-shell singularity, which is a simple pole in leading order, turns into a branch 
point as a result of the coupling of the scalar particles to an arbitrary number of soft 
photons. The propagator near the branch point exibits a behaviour which is reminiscent 
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of that of the one- loop propagator in 3+1 dimensional electrodynamics. The location of 
the branch point, which plays the role of the screening mass in the thermal problem, is 
shown to be gauge independent. However, its precise value is left undetermined by our 
present approach. 

The calculation that we have performed for SQED suggests that a similar solution 
may exist for QCD as well. If a magnetic mass exists, the mass-shell remains a pole whose 
exact location is determined by the integral equation that we have derived, assuming that 
this equation applies also to QCD; the correction to the Debye mass thus obtained is then 
identical to that calculated by Rebhan[^]. Finally, we believe that the present analysis 
should also shed light on another longstanding problem in finite temperature field theory, 
that is, the infrared singularity of the damping rates for thermal particles (see, e.g., |L{| 
and references therein). 
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Appendix 

We verify here that, close to the mass-shell k 2 — > — m 2 , the two- loop rainbow di- 
agram, Fig. 4. a, becomes as important as the one- loop graph of Fig. 2 if the infrared 
cut-off A is of the order g 2 T or less. 



A straightforward application of the Feynman rules for the action ( |3.14j ) gives (with 
(dq) = d 3 g/ (2tt) 3 ) 

V; V J J q 2 + A 2 J p 2 + A 2 [(q + k) 2 + m 2 ] 2 (p + q + k) 2 + m 2 V ; 

where we allow for a photon mass A as a convenient IR regularization, and we use the 
Feynman gauge ( = 1 for simplicity. We are interested here only in the dominant IR 
singularity of Ti^(k) as k — > ±im and A — ► 0. After working out the scalar products in 
the numerator, we isolate the terms with the largest number of factors in the denominator. 
These terms are the most singular in the double limit k 2 — > — m 2 and A — > 0. (If either 
k 2 7^ — m 2 or A 0, £( 2 ) is finite.) We denote the corresponding contribution to (|A.1|) by 
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£?i(Jfe). We have 



16m 4 (e 2 T) 2 / / _ _L _ J_ _ _ (A .2 



+ A 2 7 p 2 + A 2 [( q + k) 2 + m 2 ] 2 (p + q + k) 2 + m 2 



The terms which have been neglected in going from eq. (|A.1|) to eq. ( |A.2|) are, at most, 
logarithmically divergent in the double limit mentioned before, and do not matter for the 
power counting developed in Sect. 3. As we shall see, the integral (|A.2|) is linearly (or, 
more accurately, linearly x logarithmically) divergent in the same limit. 

To perform the momentum integrals in eq. (|A.2[) , we use the coordinate-space rep- 
resentation of the Coulomb propagators, e.g. 



k 2 + m? J Attx 
(with x = \x\) and obtain 

= / d 3 x / d 3 y e^e"™^ *— = -2a*J(k; m, A), (A.4) 



where a = e 2 T/4ir. Once the angular integrations are done, J is given by 

J(k) m, A) = ^ [°° — I" ^ sin kx e~ Xx e ^ m+x ^ Fix, y) , (A.5) 
k Jo x Jo y 

where 

F(x, y) = (l + m\x - y\) e -"*l*-»l - (l + m(x + y)) e - m{x+v) = F l + F 2 , (A.6) 

with 

^mjlx-yle-" 11 ^ 1 - (x + y) e~ m{x+y) } (A.7) 
F 2 (x,y) = e"™ 1 ^ 1 - e - m(x+y) . 

The IR divergences occuring in the original momentum integrals as k 2 — > — m 2 and A — > 
appear now as UV divergences of the integrals over x and The most singular terms are 
generated by Fi(x,y). The contribution of F2(x,y) is, at most, logarithmically divergent. 
Let Ji(k;m, A) denote the contribution of Fi to the integral (|A.5|) . For A = 0, but 
arbitrary k we have a simple expression: 

2w 2 / 4m 2 \ 
J 1 (k ] m,\ = 0) = -—^ (^TT—- 2 ~ 1 • (A.8) 
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As k 2 — > — m 2 , this has a linear x logarithmic singularity, as announced. For A 7^ 0, the 
expression of J\ is more complicated. We give here only the explicit form of the most 
singular term, valid for arbitrary k and A: 

m f m 2m + A m 2m + A 1 

J 7iJ (fc; m, A = — < —— 77 In In f . (A.9) 

«fc [m + A — ik m + 2A — ik m + A + ik m + 2A + tk) 

After inserting this in eq. ( |A.4|) , we obtain the dominant (i.e. the most singular) mass-shell 
behaviour of the two-loop diagram ( |A.1| ) as 

^ _ m 2 } ^ _ 2a 2 ™ ^ m ^ /a^ E(1)(jfe2 ^ _ m2) ^ (A1Q) 

(recall eq. ( |3.10| )). Thus, for A of the order e 2 T ~ a or smaller, and in the vicinity of 

(2) 

the mass-shell, the two-loop contribution S}^ is of the same order in a as the one-loop 
self-energy. A similar conclusion is reached in section 3.3 using power counting. 
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Figure captions 



Figure 1. Analytical structure of the screening function S(k) in the complex /c-plane, 
showing the branch point at k = im*. The contour is that used for the evaluation of the 



integral in eq. (2.2). 



Figure 2. The one- loop contribution to the self energy in the three-dimensional effective 
theory. Full line: scalar (electrostatic) field. Wavy line: transverse gluon. 



Figure 3. Approximate representation of the trajectories followed in the complex fc-plane 
by the poles of the one-loop scalar propagator. The arrows indicate the flows of the poles 
as the coupling increases. With the present choice of gauge (£ = 2), the poles become 
real when a/m^ 1.45. 



Figure 4. Two-loop contributions to the scalar self-energy which contain linear mass-shell 
divergences. Diagram (a) and (b) occur both in QCD and in SQED, while diagram (c) 
exists only in QCD. 



Figure 5. An exemple of a multi-loop contribution to the self-energy £(&), exhibiting 
power-like mass-shell divergences. 



Figure 6. The skeleton diagrams for the self-energy of the scalar propagator in SQED. 
The dominant infrared singularities are contained in diagram 6. a. 



Figure 7. Diagrammatic representation of the Dyson-Schwinger equation of SQED in the 
quenched approximation. Both the internal scalar propagator (heavy line) and the vertex 
(heavy blob) are exact quantities, while the photon propagator (wavy line) is the bare 
propagator. 
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